Cosmological dynamics of a hybrid chameleon scenario 
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We consider a hybrid scalar field which is non-minimally coupled to the matter and models 
a chameleon cosmology. By introducing an effective potential, we study the dependence of the 
effective potential's minimum and hybrid chameleon field's masses to the local matter density. In a 
dynamical system technique, we analyze the phase space of this two-field chameleon model, find its 
fixed points and study their stability. We show that the hybrid chameleon domination solution is a 
stable attractor and the universe in this setup experiences a phantom divide crossing. 
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I. 



INTRODUCTION 



Recent cosmological observations have revealed that 
our universe is currently undergoing an accelera ting 
phase of expansion P, S 1 3 |, S i i [IS EI, 
Elj El) [i3|- This transition to the accelerating phase 
has been occurred in the recent cosmological past. Sev- 
eral approaches have been proposed to explain this late- 
time acceleration of the universe. Introducing some sort 
of unknown energy component (dubbed "darfc energy"), 
with negative pressure, is one of these various approaches 
M4fU7iLi,8ilO%21^2, 23, 24, 25, 26, 27, M, M, 
[lalMEiSHelSllB. The simplest candidate for 
dark energy is the cosmological constant. However, the 
cosmological constant suffers from serious problems such 
as huge amount of fine-tuning required for its magnitude 
and other theoretical problems such as unknown origin 
and lake of dynamics [381 . Isol liol l4l| . Nevertheless, the 
standard cosmological model with a cosmological con- 
stant has no internal problems or inconsistencies both at 
the classical and quantum levels [42j. 

In this regard, the scalar fields such as quintessence 
EEHl, phantom fields [45,], tachyon fields [46., 47] and 
so on, provide a simple dynamical model for dark energy 
which can explain cosmic accelerating expansion. The 
scalar fields can directly couple to matter Lagrangian, or 
indirectly couple to the Ricci scalar |48j, |4^, l50|, |5l| . If 
there are negligible self-interactions for scalar fields, then 
the experimental bounds on such fields require it to ei- 
ther couple to the matter much more weakly than the 
gravity does, or to couple very strongly [H, [H, Is^. [55|. 
Such fields must be very light (their mass must be of the 
order of Hq, the present Hubble parameter) in order to 
evolve cosmologically today. Also, in order to have con- 
sistency with Equivalence Principal, their coupling to the 
matter must be extremely small. Chameleon cosmology 
[sgI . Wf\ is a scenario which can address this problem suit- 
ably. In the chameleon cosmology, the scalar fields evolve 
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on the cosmological time scales today, while according to 
the expectations from string theory they have couplings 
of order unity to the matter and at the same time re- 
main very light on the cosmological scales. In fact, the 
mass of the scalar field is not constant and instead, it 
depends on the local matter density. In the high den- 
sity regime, the mass of the scalar fields is large, so that 
the resulting violations of the Equivalence Principal are 
exponentially suppressed. On the cosmological scales, 
where the density is very low, the mass of the fields can 
be of the order of the present Hubble parameter (TJq) 
and cause the current acceleration of the universe. Such 
a scalar field is dubbed "chameleon" because its physical 
propert ies, such as its mass, depend on its environment 
(Sa . [STf . We note that this definition of the "chameleon" 
field is too wide and includes situations such as plasma 
density dependence of plasmon mass and also density and 
temperature dependence of elementary particles masses 
through their effective potentials. In our setup the term 
chameleon has sense only in the case of specific exponen- 
tial coupling of known quantum fields of matter to the 
chameleon field [53| . 

The present work has been organized as follows: In sec- 
tion 2, we introduce the idea of chameleon cosmology. We 
consider a quintom scalar field which is non-minimally 
coupled to the matter (a quintom field is a hybrid of a 
quintessence and a phantom field [H, [H, HI, [61I, [g^, HI] ) ■ 
In this section, by means of the conservation equation, 
we derive the matter energy density which depends on 
the hybrid scalar field. We study the effective chameleon 
potential, its minimum and the mass of two fields about 
the minimum of the potential. We show that the value of 
the scalar fields and their masses at the minimum depend 
on local energy density. In section 3, we study the cosmo- 
logical dynamics of the model in the dynamical system 
approach and we provide a detailed phase space analysis 
of the model. We find that in the parameters space of the 
model it is possible to have a chameleon dominated stable 
attractor which has a negative effective equation of state. 
We also show that in this setup, the effective equation of 
state parameter of the model crosses the phantom divide 
line and the deceleration parameter becomes negative in 
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the past. So this model has the capabihty to explain the 
late-time cosmic speed-up. 



II. 



THE SETUP 



The action of a two-field chameleon model in 4- 
dimensions can be written as follows 



J 



S = 



R 



(1) 



where = is the gravitational coupling and Cm rep- 
resents the Lagrangian density of the matter fields. Also, 
the matter fields -0''^ are coupled to the scalar fields by 
the definition g^^l = e^'^^^^i'^'^^''''' g^^, where /3(i) are di- 
mensionless constants. In this paper, we assume just a 
single matter energy density component (pm) with cou- 
pling p. Also, according to the expectations from string 
theory, we allow /? to be of the order of unity. 

In a flat Friedmann-Robertson- Walker background, 
variation of action (1) with respect to the metric, leads 
to the following equations 
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■<y5^ + -Try {4', v) + -r-Pme 



, (2) 



2H+3H 



2 2 



(3) 

To derive these equations we have assumed that the 
matter field of the universe is a perfect fluid, so that 
Pm = i^Pm- Variation of the action (1) with respect to 
the scalar fields gives the following equation of motion 

^-v3+3i?(</)-^) + ^ + ^~2At/3L.p„e«''(*+'^) =0. (4) 
dq) dip 

Equations (2)-(4) give us the energy conservation equa- 
tion of the model as follows 

Prn + 3Hp„i{l +Uj) = -K/3(l + Uj)pm{^ + 'P) ■ (5) 

The right hand side of equation (5) shows the non- 
conservation of the energy density in this setup, which 
is due to the presence of non-minimal coupling between 
the hybrid scalar field and matter Lagrangian (in other 
words, it is due to the presence of the chameleon field). 
In the absence of chameleon field (the case that /3 goes to 
zero), equation (5) simplifies to the ordinary conservation 
equation. 



I 

If we integrate equation (5), we reach the following 
expression for the matter energy density 

p„ = Ca-3(i+'^)e''^("*+'^\ (6) 
where C is a constant. This equation shows that the 
scalar fields and the matter energy density are related to 
each other via the chameleon coupling term. 

By using equations (2) and (3), we can deduce the ef- 
fective equation of state parameter in our setup as follows 



Peff ^ 02 - y-^ - 2Vi<j), if) + 2a;p^e"^('^+'^) 
Peff 4)^ - ip'^ + 2V{(t),ip) +2p„,e'^Pi't'+f) 



(7) 

Equation (4) shows that dynamics of hybrid scalar field 
in this case does not depend just only to the 1^(0, (/j), 
but instead it depends on an effective potential which is 
defined as follows 



V,f f (0, ifi) = F(0, cp) + Kfipm 6""^^*+^^ , 



(8) 



which is usually dubbed chameleon effective potential. 
This effective potential depends explicitly on the matter 
density pm ■ We note that we assume a runaway potential 
defined as 



V{(j), if) = yQe-^''("^+"^) , 



(9) 



which decreases by increment of the scalar fields. Also, 
the coupling term (^e'^Pi'f'+f)^ increases as the scalar fields 
increase (see Figure [1]). 

So, if /3 > 0, Veff has minimum. The behavior of the 
effective potential with respect to and ip is shown in 
figure [21 

We denote the values of (j> and ip, where the effective 
potential becomes minimum (the derivative of the effec- 
tive potential becomes zero), by 0mm and ipmin respec- 
tively. Then, the minimum of the effective potential in 
our model occurs at 
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FIG. 1. The behavior of the runaway potential (chromatic 
surface) and the couphng term e"'''-'^''"*'-' (the magenta meshed 
surface) with respect to the two components (</!> and if) of the 
hybrid scalar field. By increment of the scalar fields, the run- 
away potential decreases, while the coupling term increases. 



The above expression shows that, the minimum of the 
effective chameleon potential is a line. Since we deal with 
positive scalar field, ifmin leads us to a constraint on 
as follows 
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6m2a6(i+'-) 



2/3k(1 +cj) - 2^Kn- 2\/67 



(11) 



This means that if cf) satisfies the constraint equation 
(11), there is a minimum for effective potential as is 
shown in figure [5] 

The mass of the scalar field about the minimum is 
obtained from the following relation 



So, we find 



.) 



Also the mass of the scalar field tp is given by 



leading to the result 
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f Jmin 



+ «;2/32(l+^2)^^gK/3(0_„+^„„)_ 



(12) 



(13) 



(14) 



(15) 



FIG. 2. The behavior of the effective chameleon potential 
with respect to (p and ip. The effective chameleon potential 
reaches a minimum (a line in this case) during its evolution. 



As equations (13) and (15) show, (m^)mm and {m^)min 
are increasing functions of local matter density, pm ■ This 
means that, the larger values of matter density lead to 
larger values of the chameleon field's mass. 



III. 



COSMOLOGICAL DYNAMICS 



In this section, we are going to study cosmological dy- 
namics of the hybrid chameleon model introduced in pre- 
vious section in the framework of dynamical system anal- 
ysis and phase space trajectories of the model. In this 
regard, we should firstly introduce some new convenient 
dimensionless variables. These dimensionless quantities 
help us to translate our equations of the cosmological 
dynamics in the language of the autonomous dynamical 
system. In our setup, the dimensionless parameters are 
defined as follows 

K Kip K \/V 

y 



M= ^ . (16) 

By rewriting the Friedmann equation (2) in terms of the 
new variables, we reach a constraint on the parameters 
space of the model as follows 



2 I 2 I 2 

y + z +u , 



(17) 



by which we can express one of the dimensionless vari- 
ables in terms of the others. 
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By introducing a new time variable r = In a, we obtain 
the following autonomous system in our setup 



dx n n 2 ^ 2 3a; 
— — = — Sx+Smz +-— pwu + — 
or 2 2 



q o 2 V^fl 2 3a; 
dr 2 ^ 2 



— ~ —6mzx—6nzy-\ 

or ' 2 



(18) 



(19) 
, (20) 
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3u 
Y 



+ a; — y — z + lou 



.(21) 



Now, in order to analyze the cosmological evolution in 
the dynamical system approach, we should find fixed (or 
critical) points of the model. Fixed points are defined 
as those points that autonomous equations (18)-(21) all 
vanish. Eliminating z by using constraint equation (17), 
we obtain the critical points of our setup. We find four 
critical points (M, N, P, Q) and two critical lines (£1 
and £2) in our setup which we have summarized their 
properties in tables HI Hlland IIIII 

Now we discuss characters of each critical point sepa- 
rately. 

• Critical Point M: 

Point M represents either a solution with a scalar 
field's kinetic energy term domination or potential 
energy term domination, depending on the values 
of 771 and n. Also, depending on these values, this 
solution can be stable or unstable. For instance, 
by taking m = 0.6 and n = 0.5, the universe with 
this solution is an attractor, potential energy term 
dominated, meaning that if the universe reaches 
this state, it remains there forever. The value of ef- 
fective equation of state parameter corresponding 
to this solution is negative (-0.78), so this case is 
corresponding to an accelerating universe. Figure 
[3] shows the phase space trajectories of the model 
in two dimensions for /3 = 1, m = 0.6, n = 0.5 and 
w = (corresponding to potential energy domina- 
tion). The point M is shown as an attractor point 
in this plot. Also, we can see this stable point in 
3-dimensional plot as shown in figure HI If we set 
n — 0.7 and m = 1.2, the solution is a saddle, ki- 
netic energy domination. The universe during its 
evolution can reach this state but doesn't remain 
there and evolves to another state. In this case, 
the corresponding value of the effective equation of 
state parameter is positive {uJeff = 0.9) and shows 



FIG. 3. The phase space trajectories of the hybrid chameleon 
model with m = 0.6, n = 0.5, P — 1 and uj = 0. This fig- 
ure is plotted for the case with u = 0. With these values in 
the parameters space, the critical lines Ci and £2 (the solu- 
tions corresponding to the quintessence component's kinetic 
energy domination) are repeller (and so, unstable solutions). 
The critical point M (the solution with potential energy dom- 
ination) is a stable attractor. 



a decelerating universe. This saddle point of the 
phase space trajectories is shown in figure [5] This 
figure is plotted for n — 0.7 and m = 1.2. The 
point M would be corresponding to a repeller, ki- 
netic term dominated solution if for instance we set 
m = 1.6, n — 0.8 and ^ = | (see figure O. This 
unstable solution is relevant to early times cosmol- 
ogy- 

• Critical Point N: 

Critical point N exists if a; ^ 1 (in other words, 
if there is no stiff fluid in the universe). For this 
solution, the effective equation of state parameter 
is exactly the same as the equation of state param- 
eter of the matter and so the universe with this so- 
lution is effectively matter dominated. For uj < 1, 
depending on the values of m and n, this matter 
dominated solution can be stable or saddle point. 
Also, with u! — 0, this solution is always a saddle, 
independent of the values of m and n. This means 
that the universe, during its evolution, experiences 
this matter domination state and then evolves to 
another state. The point N, in both figures [Hand 
[HI is a saddle point in the phase space of the model. 
It should be noticed that since both figures have 
been plotted with oj = 0, for both m — 0.6 and 
m = 1.1, N is a saddle point. 

• Critical Point P: 

Like as the critical point N, the critical point P 
exists if 7^ 1. The properties of this critical point 
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TABLE I. Location and existence of critical points and corresponding effective equation of state parameter, m* and {xQ,yQ,UQ) 
are defined by equations (22) and (23) respectively. 



Point {x,y,u) existence 



M (m,— n,0) all m and n and all uj 

N (^rl^,#=f'l) all m and n and / 1 



p ij^, -1) all m and n and / 1 



Q (2;Q,yQ,WQ) m > m. and all a; '3 ^3 



TABLE II. Eigenvalues and dynamical characters of the fixed points. 

Point Eigenvalues ( A 1, A2, As) stability 

- fa; - vgg. ^'-v^.g"" + 3m^ - 3n^ stable if m < — (-^"/aV-^+i+^^^V-^+Dvl. 

V 2 2 2 (3n+\/6/3-3Vn2 + l){-2nfl + v/6-2flVn2 



M 



(3n+x/6/3-3yn2 + l)(-2 n|9 + v/6 -2g-^n2 + l) 

-3 + 3m2-3n^-3 + 3m2-3n2) ^ (-3n^-3-2n=^;s^-2,.^+3ny;;^)^ 

/ (3n+V6/3-3 V"^ + l)(-2"/9+^-2/3V"-^ + l) 

|- -3+3^^+2^v/6/3^^-2nV6ga, ^ ^^^^^j^ if < 1 and m < (3-3^ V2nv^^^) 



+ fo), - | + ftj^ saddle point if = 



-3+3^V2,n^/3^-2n^/3^ Stable ff < 1 and m < (3-3^V2nv^g^)v^ 



2"^' ~l + I'- 



-| + |w, — I + fcj ) saddle point if cj = 



Q (-^iQ, A2Q, Asq) stable 




FIG. 4. 3-dimensional phase space trajectories of the model 
with m — 0.6, n — 0.5, P ~ 1 and uj = 0. For these values 
of parameters, we have a potential energy dominated solution 
(point M) which is a stable attractor. The points N and P 
(corresponding to the effectively matter dominated era) are 
saddle points. 



FIG. 5. 3-dimensional phase space trajectories of the model 
with m = 1.2, n = 0.7, P — 1 and lj — 0. For these values of 
parameters, the potential energy dominated solution (point 
M) is a saddle point, while the point Q (corresponding to the 
chameleon dominated solution) is a stable attractor solution. 



figure has been plotted for positive u while P is 
located at it = —1. 



is the same as the properties 01 the critical point 
N. In figure SI we can see the point P as a saddle 
point in the phase space trajectories of the model. 
This point is not shown in figure |31 because this 



• Critical Point Q: 

The critical point Q exists if there is the following 
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constraint on the parameters space of the model 



FIG. 6. The phase space trajectories of the model with m = 
1.6, n = 0.8, /3 = 1 and uj — ^. Like as figure 3, this figure 
is plotted also for the case with u = 0. With these values 
of parameters, in contrast to figure 3, the critical line C\ is 
stable attractor, while the critical line C2 remains unstable. 
Also, in this case the critical point M is a repeller. 



^6/32^2 + 72 + 144n2 + 72w - 24V6/3wn 
m > —f3uj + — = TO* . (22) 

I 

If this critical point exists, its location is at 

{xQ, yq, uq) (23) 

where 

I 



_ 127i2^w - 3m^/6 + 6a;/3 - 12mn;3w + im^/QuJ'^ + Gw^/J _ 2^/QP'^uJ'^n + 2^6/32^2^ 
-\/6(— Stow + 6to — 6naj + 6n + -\/6a;2/3 — 3-y/6i^/3)(?T. — to) 

_ + 6/3a;2 - 2%/6/32a;2n - 3%/6n + 3%/6nw2 _^ I2ml3ujn + 2V6/32a;2^ _ 12to2w;3 
^'^ \/6(6?T. + 6to — Qnu} — Qmuj — 3/3aj\/6 + j3ijj'^\/&){ra — n) 

and 

\/3ri + 3to — 2/3a;\/6 — Qnuj — 6muj + 2\/6a;2/3 + 3nw2 + Suj'^m — 2(3'^uj'^n + 
(—6 too; + 6 to — 6 nuj + Qn + \f^uj^fi — 3 P u V6) {n — to)^ 



X y -12to2 + 6 + 12n2 + 6w + 2\/6tow/3 - 2V6ujl3n. 

I 

Also, (AiQ, A2Q, Asq) are the eigenvalues of the fol- lowing Matrix 
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I dx 



M = 




(24) 



{x.,y.u) = {xQ.,yQ,UQ) 



where prime refers to derivative with respect to r 
(see equations (18)-(21)). We note that since the 
eigenvalues of this solution are so lengthy and com- 
plicate, we avoid to express them here explicitly. 
The universe with this solution experiences an ac- 
celerating phase if 



and n = 0.7. With these values of the parameters, 
the critical points M and N are saddle points and 
the points Q is a stable attractor. In summary, 
we can say that with these parameters values, the 
universe during its evolution reaches a kinetic en- 
ergy dominated era and then evolves to a matter 
dominated era. After that the universe evolves to 
a chameleon dominated era and remains there for- 
ever. Also, with this choice of values, the value of 
the effective equation of state parameter is about 
— 1.1. This means that the universe in this param- 
eters space experiences late time acceleration and 
its stable state lies in phantom-like phase. 




)(W+ -) < XQ 



< 



^-ul{u: + -). (25) 



• Critical Line ci'- 

In this model, we also have two critical lines. 
The critical line £i, which is located at 
(a:*, yjxl — 1, 0), is a quintessence's kinetic term 
dominated solution. So, there is a constraint on 
this solution. This kinetic energy dominated solu- 
tion exists if 



Point Q of figure [5] is a stable chameleon domi- 
nated solution and this means that if the universe 
reaches this state, remains there forever. In this 
solution, the universe can experience the late time 
acceleration. Figure [5] shows 3-dimensional phase 
space of our setup with /?=l,a; = 0, m = 1.2 



> 1 or a;, < 1 . (26) 

The stability of this line depends on the model pa- 
rameters values. In the parameters space which 
satisfies the following constraint, this solution is a 
stable attractor 



TO > 



\/6n(3 CO x'^ + VGxV—^ + x'^n(3 ui — /3uj \/6x — \/6n(3 lu — fiuj \/6\/— 1 + x'^ 
x{l3uj VQx + /3 cj VQV-I + x^ + 3 cj - 3) 



-1-3 -I- 3 nV— 1 + x'^u} — 3 7i\/— 1 + x'^ — 3 w 
x{Puj V6x + /3uj VeV-l + + 3 w - 3) 



(27) 



In figure |3] where the chosen values of the parame- 
ters do not satisfy the constraint equation (27), the 
critical line Ci is unstable. But in figure |6]this line 
is a stable attractor because the values of parame- 
ters used to plot this figure satisfy the mentioned 
constraint. 

• Critical Line £2- 

This phantom's kinetic energy dominated critical 
line is located at (— x*, —\/x1 — 1, 0). So, for ex- 
istence of this solution the constraint equation (27) 
should be satisfied. In contrast to £i case, the crit- 
ical line £2 is an unstable solution for all values of 
the model parameters space (see figures [3] and IH]). 

A parameter which gives us a suitable background 



to understand the dynamics of the universe, the na- 
ture of dark energy and the possibility of crossing of the 
phantom-divide line, is the equation of state parameter. 
According to the recent observational data, the equation 
of state parameter of the dark energy crossed the phan- 
tom divide line (w = —1) in the near past. It is shown 
that, considering a dynamical dark energy component in 
a cosmological setup enables the model to explain these 
observational evidence 



|16L 17 IS 




|2J, |25|, |26l 1271 128|, |29j, ISO, |31|, 13 
[7]shows the evolution of We/ / with respect to the red-shift 
parameter in our setup. This figure has been plotted by 



adopting the ansatz with a 

-at 



lit 



ipo e~ 



and 

pQe^'"- {a, 5 and j5 are positive constants). Note 
that in this figure, the values of parameters are the same 
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TABLE III. Location, eigenvalues and dynamical characters of the critical lines. Note that, mci is defined in (27). 



Line 


{x,y,u) 




existence 


Eigenvalues 




stability 


£1 


(Xtj^/xl — 


1,0) 


all m, n and a; if 2; > I 


(AlCi,A2£i, ^3Ci) 


I 


stable if 
m > mci 


C2 


{x,,-y^xi - 


-1,0) 


all m, n and a; if x < — 1 




1 


unstable 





FIG. 7. The evolution of the effective equation of state pa- 
rameter of the hybrid chameleon model versus the red-shift 
for the case that m = 1.2, n = 0.7, /3 = 1 and to = 0. 



FIG. 8. The evolution of the deceleration parameter of the 
hybrid chameleon model versus the red-shift for the case that 
m = 1.2, n = 0.7, /3 = 1 and a; = 0. 



as figure [5] (where the chameleon dominated solution is 
a stable attractor). Also, we have set ap — I, <j)o — 0.8, 
a = 1, 6 = 1.6 and (3 — 1. As figure shows, in the 
presence of the chameleon field, the universe enters the 
phantom phase in the near past at z ~ 0.26. So, in this 
model, the universe experiences a smooth crossing of the 
phantom divide, We// = — 1, line. 

The deceleration parameter, q, is another important 
parameter in cosmological evolution. Its positive value 
(corresponding to a < 0), shows the decelerating expan- 
sion of the universe and its negative value (corresponding 
to a > 0), shows the accelerating expansion, q is defined 
as follows 



H 



1. 



(28) 



The evolution of the deceleration parameter versus red- 
shift is plotted in figure El We can see from this figure 
that the deceleration parameter has become negative in 
the past at z ~ 0.84, meaning that the universe has en- 
tered to an accelerating phase at z ~ 0.84. So, a hybrid 
chameleon model has the phantom like behavior and can 
explain the late time cosmic acceleration of the universe 
in an observationally viable manner. 



IV. SUMMARY 

In the current paper, we have studied a hybrid 
chameleon model in details. By choosing a runaway po- 
tential and by using of the coupling term e'^^^'^^^\ we 
have introduced an effective chameleon potential. In con- 
trast with runaway potential and coupling term, this ef- 
fective chameleon potential has a minimum. We have 
derived this minimum and the messes of the hybrid 
chameleon field in the minimum. We have found that 
the minimum of the effective chameleon potential is a 
line and in order to existence this minimum, there is a 
constraint on the parameters space of the model. It has 
been shown that the chameleon field and the matter en- 
ergy density is related to each other. Also, we have shown 
that the larger values of the density leads to larger values 
of the chameleon field's masses. Then we have consid- 
ered the cosmological dynamics of the hybrid chameleon 
model. In a dynamical system approach, we have studied 
the phase space trajectories of the model and its stabil- 
ity. We have found four critical points and two critical 
lines in our setup and have derived their eigenvalues in 
order to find the stable solutions. The critical lines and 
points can be attractor, repeller and saddle, depending 
on the values of the parameter spaces. By analyzing two 
and three dimensional phase spaces, it was shown that 
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there are matter domination, scalar field's kinetic term 
domination and chameleon domination stable solutions, 
depending on the values of the parameters space. With 
some values of the parameters space, we have found the 
solution where the universe after passing the scalar field's 
kinetic term dominated era (repeller) and the potential 
energy dominated and matter dominated era (saddle), 
reach a stable attractor chameleon dominated era and 
remains there forever (one can see the saddle point and 



attractor in figure [5]). Also, the effective equation of 
state parameter in this case is below —1. By adopting 
the parameter's values used in figure [5l we have ana- 
lyzed the late time dynamics of the hybrid chameleon 
model. We have shown that the universe in this setup 
has crossed the phantom divide line at z ~ 0.26 and has 
entered to the accelerating phase at z ~ 0.84. So, a hy- 
brid chameleon model can explain the late time cosmic 
acceleration and realize the phantom behavior. 
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